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ABSTRACT:  The  problem  of  vibrational  relaxation  of  nitrogen  in 
Hypersonic  Tunnel  No.  4  at  the  U.  S.  Naval  Ordnance  Laboratory 
has  been  analyzed  for  supply  conditions  of  3600°R,  3960°R,  5400°R, 
and  pressures  of  80  and  100  atmospheres.  A  comparison  between 
equilibrium,  relaxing,  and  frozen  flows  has  been  made  for  flow 
in  the  nozzle  and  around  a  pitot  probe  placed  in  the  flow.  The 
results  are  shown  as  plots  of  static  pressure,  static  temperature, 
free  stream  velocity,  and  stagnation  pressure  ratio  versus 
nozzle  area  ratio.  It  has  been  found  that  the  flow  freezes  at 
very  short  distances  downstream  of  the  throat.  When  plotted 
versus  nozzle  area  ratio,  the  ratio  of  impact  to  stagnation 
pressure  gives  identical  results  in  the  relaxing  and  frozen  flow 
cases.  The  authors  recosmend  that  frozen  flow  be  used  as  a 
convenient  and  reasonably  accurate  assumption  for  determining  the 
test  section  conditions  in  NOL  Hypersonic  Tunnel  No.  4. 
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This  is  a  report  of  a  study  under  takes  to  compare  the  equilibrium , 
relaxing,  and  frozen  flow  in  a  conical  nozzle  of  Hypersonic 
Tunnel  No.  4  at  the  U.  S.  Naval  Ordnance  Laboratory.  Of  partic¬ 
ular  Interest  was  the  interpretation  of  the  pressure  recorded  by 
a  pitot  probe  placed  in  the  flow. 
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SYMBOLS 

area  of  nozzle  section,  ft^ 
speed  of  sound,  ft/sec 
defined  in  equation  (10) 
theraal  speed  of  aolecule,  ft/sec 
average  speed  of  a  aolecule,  ft/sec 
specific  heat  at  constant  pressure 
internal  energy  of  a  gas 
enthalpy 

Planck's  constant 
Boltzaann's  constant 
aean  free  path,  ft 
Mach  miaber 

nuaber  of  collisions  required  for  equllibriua 
pressure,  ataospheres 

universal  gas  constant,  1776  ft/lb/slug°R 
entropy 

teaperature,  °R 
tlae,  sec 
velocity,  ft/sec 

distance  froa  nozzle  throat  downstreaa,  ft  or  in  as 
noted 

standoff  distance  between  shock  and  tip  of  pitot  probe 
ratio  of  specific  heats 
coefficient  of  viscosity,  slug/ft/sec 
vibrational  frequency 
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0T  characteristic  vibrational  tesperature,  °R 

p  density,  lb/ft^ 

t  relaxation  tise,  sec 

Subscripts 

1  indicates  conditions  before  shock  wave 

2  indicates  conditions  after  shock  wave 

F  pertains  to  frozen  flow 

o  stagnation  or  total 

r  pertains  to  rotational  sode 

t  pertains  to  translational  node 

v  pertains  to  vibrational  sode 

*  throat  conditions 
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INTRODUCTION 

Hypersonic  Tunnel  No.  4  at  the  U.  S.  Naval  Ordnance 
Laboratory  has  recently  been  Modified  by  the  addition  of  a 
graphite  heater.  This  tunnel  operates  at  a  typical  supply  con¬ 
dition  of  4000°R  and  100  atmospheres,  achieving  Mach  nusbers 
above  15.  It  uses  nitrogen  rather  than  air  as  the  test  gas. 

At  the  above  supply  condition  a  considerable  amount  of  the 
Internal  energy  of  the  gas  is  found  in  the  vibrational  energy 
state.  Figure  1  shove  the  percentage  of  the  internal  energy 
found  in  the  vibrational  state  as  a  function  of  temperature. 

It  also  shovs  the  vibrational  energy  as  a  percentage  of  the 
enthalpy.  These  curves  are  obtained  from  equations  to  be  given 
later  in  this  report. 

Unlike  the  translational  and  rotational  energy  states,  this 
vibrational  energy  requires  a  considerable  number  of  lntermolec- 
ular  collisions  for  adjustment  to  any  new  state.  In  this  process 
of  adjustment,  one  can  consider  tvo  limiting  cases:  equilibrium 
flow  and  frozen  flow.  In  the  former  the  energy  in  the  vibrational 
degrees  of  freedom  depends  only  on  the  translational  temperature. 
In  the  latter,  the  energy  is  locked  in  the  vibrational  state  and 
does  not  change. 

The  purpose  of  this  report  is  to  consider  some  effects  of 
vibrational  energy  on  the  flow  in  our  hypersonic  test  facility. 

We  have  analyzed  the  tvo  limiting  cases  outlined  above.  We  have 
also  considered  the  gas  to  be  relaxing,  that  is,  adjusting  its 
vibrational  energy  at  a  finite  rate.  In  the  first  part  the 
nozzle  flov  is  considered.  The  second  part  discusses  the  flov 
around  a  pitot  probe  placed  in  the  flov,  and  finally  the  results 
of  the  first  tvo  parts  are  combined  to  give  useful  information  on 
the  interpretation  of  experimental  pitot  pressure  results  in 
Tunnel  No.  4. 

We  have  neglected  ionization  and  dissociation  since  these 
are  not  expected  for  our  flov  conditions. 


NOZZLE  FLOW 


A.  Equilibrium  Flov 

The  internal  energy  of  a  gas  is  composed  of  a  translational, 
rotational,  and  vibrational  mode: 

E  ■  +  Er  +  Ev 
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or,  combining  translational  and  rotational: 

E  -  |  RT  ♦  Ev  (1) 

The  vibrational  energy  (Ev)  associated  with  the  diatomic 
molecule  Is  usually  considered  as  that  of  a  simple  harmonic 
oscillator.  From  reference  (1),  when  the  flow  is  In  complete 
thermodynamic  equilibrium,  Ey  may  be  expressed  as  a  function  of 
temperature  in  the  following  way: 

E’ ' R  JW&tl  <2> 

where  0y  is  a  characteristic  temperature  for  a  particular  gas 
which  depends  on  Planck's  constant  (h),  Boltzmann's  constant  (k) , 
and  vibrational  frequency  (v)  in  the  following  way: 


For  diatomic  nitrogen,  6y  "  6102°R. 

The  equilibrium  flow  in  the  nozzle  is  treated  assuming  isen 
tropic  flow.  This  is  equivalent  to  using  the  usual  one- 
dimensional  Euler  equations,  as  seen  in  the  following  way: 

Consider  the  energy  equation: 

H  +  ^  ua  -  constant  (3) 


where  H  -  E  ♦  P/p. 

Differentiation  yields: 

dH  -  -  u  du 

On  the  other  hand,  the  definition  of  entropy  yields 


(4) 


dS  -  -yCdE  +  P  d  i]  -  ^CdH  -  i  dP] 


(5) 


For  isentropic  flow,  dS  -  0.  Therefore: 

dH  -  i  dP 
P 

Using  equation  (4)  one  finally  obtains: 

u  du  -  -  -  qP 

P 


(6) 


2 


NOLTR  63-221 


which  is  the  usual  equation  for  one-dimensional  inviscid  flow. 
Thus,  the  two  methods  are  equivalent. 

Let  us  proceed  now  to  obtain  the  entropy  equation.  Consider 
first  the  equation  of  state: 


P  -  pRT  (7) 

In  reality.  P  ■  pRTZ,  where  Z  is  the  compressibility  factor  and 
depends  on  the  virial  coefficients.  However,  from  the  data  in 
reference  (2) ,  a  plot  was  made  of  P  versus  T  for  constant  values 
of  Z  (fig.  2)  showing  that  for  pressures  less  than  100  atmospheres, 
the  error  in  assuming  Z  -  1  is  very  small.  We  have  thus  set 
Z  -  1  and  used  equation  (7). 

By  means  of  equations  (1).  (3)  and  (7)  one  obtains: 

H  ■  |  RT  +  Ev  (3) 

Note  that  H  is  a  function  of  T  only.  From  equations  (5)  and 
(3): 

dS  1  dH  dP 

TITHT'T 


Integration  yields: 

§  "  J  KT  HT  dr  -  lnP  +  (constant  of  integration) 


Substituting  (8)  into  this  equation  and  integrating,  one  obtains: 

^  ■  J[ln  T  +  1]  -  In  P  -  ln[l  -  e“6v^r]  +  t[  §v/T — ]  +  constant 

where  P  is  in  atmospheres  and  T  in  degrees  Rankine.  The  value 
of  the  constant  was  obtained  by  matching  this  formula  to  the 
NBS  data  (ref.  (3))  at  one  atmosphere  and  540°R.  We  obtained 
-2.4767  for  pressure  in  atmospheres  and  temperature  in  degrees 
Rankine.  Thus,  the  entropy  equation  is: 

|  -  |[ln  T  +  1]  -  In  P  -  ln[l  -  e"e»/T]  +  -f[ ]  *  2-«?67 

e  —  1 

(9) 
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The  Bpted  of  sound  in  equilibrium  flow  is  now  found  in  a 
straightforward  but  tedious  Manner: 

*’  ■  (sf)s  * 


Differentiation  of  (7)  and  (9)  according  to  the  expression  above 
gives  for  the  speed  of  sound: 


where 


The  equilibrium  flow  in  the  nozzle  is  then  found  as  follows: 
The  entropy  equation  (9)  is  used  to  find  S  in  the  stagnation 
chamber.  The  static  pressure  in  the  nozzle  as  a  function  of 
temperature  is  found  from  the  same  equation  using  this  constant 
value  of  S.  The  throat  conditions  are  found  by  maximizing  the 
Mass  flow  per  unit  area,  pU.  At  any  station  along  the  nozzle, 
the  area  ratio  A/A*  is  found  from  the  conservation  of  mass 
equation: 


The  results  of  this  section  for  a  typical  supply  condition 
are  shown  in  figures  3,  4,  and  5  labeled  "equilibrium."  These 
figures  show  pressure,  temperature  and  velocity  as  a  function  of 
the  area  ratio  of  the  nozzle.  The  other  two  curves  appearing  on 
these  figures  will  be  discussed  later. 


Partial  Equilibrium 


Vhereas  in  the  case  of  couplets  equilibrium  tbs  vibrational 
energy  depends  only  on  ths  translational  temperature  (eq.  (2)), 
in  ths  partial  equilibrium  cans  this  squatlon  no  longsr  appliss. 
Instsad,  there  exists  an  additional  differential  equation  (the 
rat#  equation)  for  determining  1T. 
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In  obtaining  the  rate  equation,  it  is  convenient  to  define 
two  new  variables:  the  relaxation  time,  t,  and  the  vibrational 
temperature,  Tv.  The  concept  of  relaxation  time  is  generally 
employed  to  characterize  the  time  required  for  the  adjustment  of 
the  various  degrees  of  freedom  to  equilibrium.  It  represents 
the  time  required  to  approach  1  -  1/e  of  the  final  equilibrium 
value.  From  reference  (4)  the  relaxation  time  in  seconds  as  a 
function  of  temperature  and  pressure  is: 

✓a 


T 


“0  •  46  l  T 

- nr 


(id 


where  T  is  in  degrees  Rankine  and  P  is  in  atmospheres. 

Equation  (11)  is  shown  graphically  in  figure  6.  It  is 
noted  from  this  figure  that  the  relaxation  time  increases  as  the 
gas  travels  through  the  nozzle,  since  both  pressure  and  tem¬ 
perature  are  decreasing.  As  a  result,  it  is  found  in  hypersonic 
nozzles  that  the  vibrational  energy  is  unable  to  follow  these 
rapid  changes.  Therefore,  the  vibrational  energy  can  no  longer 
be  obtained  from  equation  (2).  However,  one  can  always  find  a 
temperature  such  that  this  equation  be  true.  This  temperature 
is  known  as  the  vibrational  temperature,  Tv.  Thus,  the  vibra¬ 
tional  energy  in  the  case  of  partial  equilibrium  is: 

Ey  *  R  •  a  (12) 

e °v/Ty_1 

The  rate  equation  may  now  be  introduced.  From  reference  (5) 
this  is: 

^  ■  ^[Ey(T)  -  Ey(Iy)"l  (13) 

It  will  be  noted  that  the  right-hand  side  is  proportional  to  the 
difference  between  two  vibrational  energies — the  actual  energy 
and  the  energy  the  gas  would  have  if  it  were  in  equilibrium  with 
the  translational  energy. 

The  flow  through  the  nozzle  may  now  be  solved  by  adding  the 
rate  equation  (13)  to  the  equations  used  in  the  equilibrium  case 
and  manipulating  them  to  obtain  three  differential  equations  in 
u,  r,  and  Tv  which  are  solved  simultaneously. 

Summarizing,  we  have  the  following  equations: 

Momentum:  u  -  -  1  ||  (6) 
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Continuity:  puA  ■  constant 

Energy:  H  +  %  ua  -  constant 

State:  P  -  pRT  (7) 

Enthalpy:  H  ■  ^  RT  +  Ev  (8) 

Vibrational  Energy  0 

for  Equilibrium:  Ev  -  R  -  n  V™ —  (2) 

e0v/T-l 

Actual  Vibration  Energy:  Ev  -  R  — -  (12) 

eev/rv.i 

Rate  Equation:  ^  -  1[ev(T)  -  EV(TV)]  (14) 

where  t  is  given  by  (11). 

The  foregoing  system  of  equations  which  include  the 
vibrational  rate  equation  can  be  reduced  to  a  set  of  three  dif¬ 
ferential  equations  each  of  which  contain  only  the  variables  T,  Tv, 
u,  and  x  for  a  given  nozzle  with  a  fixed  stagnation  condition. 
These  three  equations  are: 


35  "  Vu  A  RT  [x  a;  +  75t(Ev<t)  ■  Ev(TV>)]  (15) 
T  u 


dTv  _  1  EV(T)  -  EV(TV) 

air  “  uEv'(Tv)  - t 

(16) 

dr  .  -ar  du  _  ...  .  dTyi 

Hx  7Rlu  Hx  +  Ev  (Tv)  3x  J 

(17) 

where  the  prime  Indicates  differentiation  with  respect  to  Tv. 
in  this  report  we  shall  use  an  equation  for  the  nozzle  area 
which  approximates  that  used  in  Hypersonic  Tunnel  No.  4.  It  is: 

-  1  +  700xa 

(A  simple  calculation  from  this  formula  shows  the  radius  of 
curvature  at  the  throat  to  be  3.72  in.  This  is  the  value  used 
in  the  actual  nossle.) 
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We  haw  therefore: 


1  dA  1400x 

j  Hx  "  r;  700x3 

where  x  is  any  distance  downstream  of  the  throat  in  feet. 

These  equations  were  solved  simultaneously  on  the  IBM  7090 
computer  at  NOL.  It  was  assumed  that  the  vibrational  energy  is 
in  equilibrium  with  the  translational  and  rotational  energy 
modes  for  the  flow  process  from  the  nozzle  stagnation  conditions 
up  to  the  throat.  Downstream  of  the  throat  the  vibrational 
energy  is  allowed  to  relax. 

That  this  assumption  is  reasonable  is  seen  in  the  following 
way.  If  t  is  the  average  time  the  gas  takes  in  traversing  the 
upstream  region,  the  ratio  t/T  gives  us  an  indication  of  whether 
the  gas  is  in  equilibrium  or  not.  For  example,  if  t/T  »  1,  the 
relaxation  time  will  have  elapsed  before  the  gas  arrives  at  the 
throat,  and  thus  the  flow  will  be  in  equilibrium.  With  this  in 
mind,  we  have  plotted  a  curve  of  t/T  versus  T  for  values  of 
constant  pressures  (fig.  7).  For  t  we  have  used  t  -  x/u,  where 
x  is  the  distance  from  the  start  of  the  flow  to  the  throat  (.75 
in.)  and  u  -  a*/2  where  a*  is  the  speed  of  sound  at  the  throat, 
a  function  of  the  stagnation  conditions. 

The  figure  shows  that  for  the  range  of  pressures  greater 
than  forty  atmospheres,  t  is  greater  than  10t  for  temperatures 
above  3400°R.  Thus,  the  assumption  of  equilibrium  flow  up  to 
the  throat  is  indeed  reasonable  for  the  conditions  in  Hypersonic 
Tunnel  No.  4. 

The  results  of  the  partial  equilibrium  solution  are  shown 
in  figures  3,  4,  and  5  as  curves  of  static  pressure,  static  tem¬ 
perature,  and  free  stream  velocity  versus  area  ratio.  In 
figure  4,  we  have  also  plotted  Tv,  and  the  curve  indicates  that 
freezing  of  the  vibrational  energy  occurs  at  a  short  distance 
from  the  throat.  For  the  case  shown,  freezing  occurred  at 
x  -  0.05  ft.  The  freezing  point  has  been  chosen  as  that  point 
where  the  vibrational  temperature  has  completed  95  percent  of  its 
total  change.  The  total  nozzle  length  used  in  Hypersonic  Tunnel 
No.  4  is  3  ft.  from  throat  to  exit.  For  other  stagnation  condi¬ 
tions  up  to  5400°R  and  100  atmospheres,  freezing  occurred  at 
values  of  x  less  than  0.055  ft.  In  figure  3  we  have  plotted  the 
final  frozen  vibrational  temperature  (Tyj.)  as  a  percentage  of 
the  supply  temperature  for  various  stagnation  conditions.  It  is 
seen  that  64%  <  Typ/T0  <  77%  for  40  atm  <  P  <  100  atm  for 
3600°R  <  V  <  5400°K.  This  ratio  decreases  as  pressure  and  tem¬ 
perature  increase. 


7 


NOLTR  63-221 


Thus  we  see  that  for  the  most  part  of  this  3  ft.  nozzle, 
the  vibrational  energy  is  completely  frozen.  We  shall  consider 
next  the  completely  frozen  flow  case. 

C.  Frozen  Flow 


The  other  limiting  case  in  high- temperature  flow  is  the 
case  in  which  the  vibrational  energy  is  constant  throughout  the 
nozzle.  This  is  known  as  frozen  flow.  Since  the  vibrational 
energy  does  not  participate  in  the  flow  process,  this  case  can 
be  considered  identical  to  perfect  gas  flow.  Mathematically, 
one  can  combine  the  one-dimensional  flow  equation  with  the 
expressions  for  enthalpy  and  manipulate  them  to  obtain  P/pV 
constant.  This  shows  frozen  flow  to  be  equivalent  to  perfect 
gas  flow.  Reference  (6)  gives  a  mathematical  proof  of  this  fact. 

We  have  used  the  same  procedure  in  solving  this  flow  as  was 
used  for  the  equilibrium  case,  this  time  making  dGv  “  0.  The 
expression  for  entropy  is  now: 

|  -  3.5  In  T  -  In  P  -  2.4767 

The  results  are  shown  in  figures  3,  4,  and  5  as  curves  of  static 
pressure,  static  temperature,  and  free  stream  velocity  versus 
nozzle  area  ratio. 


PROBE  FLOW 


A .  Shock  Wave 

Vibrational  relaxation  in  the  region  near  a  pitot  tube  may 
affect  the  pressure  readings  recorded  by  the  probe.  The  tem¬ 
perature  behind  the  shock  wave  created  by  the  probe  increases 
as  Mxa  and  becomes  large  enough  so  that  the  molecular  collisions 
excite  the  vibrational  degrees  of  freedom  of  the  diatomic  mole¬ 
cules.  Immediately  behind  the  shock  the  vibrational  relaxation 
process  will  start.  In  order  to  determine  the  effects  of  this 
process  on  the  pressure  readings,  we  will  consider  the  flow  as 
it  passes  through  the  shock  and  across  the  shock  layer. 

As  the  gas  moves  through  the  shock  wave,  no  vibrational 
relaxation  is  expected  (the  flow  is  frozen)  since  the  number  of 
molecular  collisions  within  the  shock  wave  is  small  compared  to 
the  number  required  for  equilibrium.  This  may  be  shown  as 
follows:  An  approximate  figure  for  the  number  of  collisions 
required  for  equilibrium  can  be  obtained  if  we  know  the  thermal 
speed  of  the  molecule  (C),  the  mean  free  path  (L)  and  the 
relaxation  time.  From  reference  (7): 
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c  -  yzirr 


L 


rhe  viscosity 


of  nitrogen  is  found  from  Sutherland's  relation: 


2.14  x  10*# 


T*> 

TV  150 


ft-sec 


with  T  in  °K.  Since  the  relaxation  time  is  known  (eq.  (11)) , 
the  number  of  collisions  required  for  the  excited  gas  to  come 
to  equilibrium  is: 


N 


1.098  x  lO^Cf  +  180]  e‘°,4tl 

13  T3/* 


(18) 


It  is  noted  that  N  depends  only  on  temperature.  This  equation 
has  been  plotted  in  figure  (9)  where  it  is  seen  that  N  is  very 
large,  from  30,000  at  T  -  5000°R  to  700,000  at  r  -  2000°R.  On 
the  other  hand,  reference  (3)  tells  us  that  the  transition 
across  a  shock  wave  from  temperature  Ti  to  Ta  takes  place  over  a 
distance  of  the  order  of  one  mean  free  path.  That  is,  a  mole¬ 
cule  undergoes  only  a  few  collisions  in  its  passage  through  a 
shock  wave,  and  thus  there  are  not  enough  of  them  for  any 
relaxation  to  occur. 

B.  Shock  Layer 


As  the  gas  travels  through  the  shock  layer,  the  relaxation 
process  begins.  However,  because  of  the  very  short  standoff 
distance  between  the  shock  and  the  probe,  the  gas  can  be  assumed 
to  remain  frozen.  An  indication  of  this  fact  may  be  obtained 
by  calculating  t/r,  the  ratio  of  transit  time  to  relaxation 
time,  for  this  region.  If  this  ratio  is  very  small,  the  flow 
can  be  assumed  to  be  frozen.  The  transit  time  for  the  shock 
layer  is  given  by 


t 


(19) 


where  A  is  the  standoff  distance.  From  reference  (9),  for  a 
strong  shock  wave  where  pa/px  -  6,  the  standoff  distance  is 
approximately  0.0011  ft.  when  the  tip  of  the  pitot  probe  is 
approximated  by  a  sphere  of  radius  equal  to  1/8  in.  The  value 
of  ua  is  an  average  value  of  the  velocity  behind  the  shock.  The 

relaxation  time  is  obtained  from  equation  (11)  using  the 
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temperature  and  pressure  behind  the  shock,  Pa  and  Ta .  The 
variables  u3 ,  Pa ,  and  Ta  are  obtained  from  the  supply  conditions 

assuming  completely  frozen  flow  from  reservoir  to  shock.  This 
Is  permissible  since  we  are  Interested  In  obtaining  only  the 
order  of  magnitude  of  the  ratio  t/t.  Hence,  given  the  supply 
conditions,  we  can  use  perfect  gas  tables  such  as  reference  (10) 
and  obtain  ua ,  Pa ,  and  Ta  and  thus  t,  t,  and  t/T.  This  has  been 
done  for  a  typical  operating  Mach  number  in  Hypersonic  Tunnel 
No.  4  (M-17)  and  in  figure  10  we  have  plotted  the  ratio  t/r  in 
the  shock  layer  for  various  supply  conditions.  It  is  seen  that 
for  the  range  3S00°R  <  T  <  5000°R  and  40  atm  <  P0  <  100  atm. 

The  ratio  t/T  is  in  the  range  of  10“4  to  10“3.  Thus,  we  expect 
frozen  flow  in  the  shock  layer. 


COMPLETE  FLO* 

The  three  cases  previously  discussed  separately  are  now 
combined  and  compared  using  the  ratio  of  the  supply  stagnation 
pressure  (P0i)  to  the  pressure  read  by  the  pitot  probe  (P0a). 
For  conditions  before  the  shock  we  have  used  the  nozzle  flow 
data  leading  up  to  figures  3,  4,  and  5  and,  hence,  have  assumed 
the  nozzle  flow  to  be  in  equilibrium,  partial  equilibrium,  and 
frozen.  For  the  flow  through  the  shock  and  the  shock  layer,  we 
have  assumed  the  two  limiting  cases:  equilibrium  and  frozen. 
The  final  result  for  each  case  is  a  curve  showing  the  variation 
of  Poa/Poi  with  Mach  number  and  area  ratio  for  various  typical 
stagnation  conditions.  The  cases  are  outlined  in  TABLE  1. 


TABLE  1 


CASE 


UPSTREAM  DOWNSTREAM 

OF  THROAT  OF  THROAT 


SHOCKWAVE  SHOCK 

LAYER 


1 

2 

3 


Frozen 

Equilibrium 

Equilibrium 


Frozen 

Relaxing 

Equilibrium 


Frozen 

Frozen 

Equilibrium 


Frozen 

Frozen 

Equilibrium 


CASE  1:  Since  frozen  flow  is  identical  to  perfect  gas  flow,  Cp 
is  a  constant  throughout,  and  the  stagnation  pressure  ratio  may 
be  obtained  by  means  of  the  NACA  Tables,  reference  (10). 


CASS  2:  Co  is  still  constant  across  the  shock  and  shock  layer, 
and  reference  (10)  may  also  be  used.  However,  the  conditions  in 
region  1  (immediately  before  the  shock)  must  be  obtained  from 
the  partial  equilibrium  nozzle  flow  data. 
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CASF  3:  Cp  is  now  not  constant,  and  reference  (10)  may  not  be 
used .  Instead,  we  shall  write: 


Poa  .  Pi  Pa  Ppa 

Poi  Poi  Pi"  Pa 

Pi 

a.  The  ratio  ■*— is  obtained  from  the  equilibrium 

P0i 

flow  data  leading  to  figure  3. 


(20) 

nozzle 


Pa 

b.  The  ratio 
outlined  below:  1 


is  found  using  an 


iteration  procedure 


According  to  reference  (11),  the  velocity  ratio  -p  is  a 

good  indicator  of  the  equilibrium  state  of  the  gas.  1  The 

value  of  this  ratio,  it  is  shown,  decreases  as  additional 
degrees  of  freedom  are  excited.  The  shock  equations  may  be 
expressed  in  terms  of  this  ratio,  and  once  its  value  is  found 

(by  iteration  on  tht.se  equations)  they  may  be  solved  to  obtain 

the  values  of  the  variables  behind  the  shock.  The  iteration 
procedure  is  as  follows: 


The  conservation  equations  across  the  shock  are: 


P 1U1  -  Pa  Ua 

Pi  +  Pi“ia  -  Pa  +  Pa  Ua* 
Hi  +  Ux 3  •  Ha  +  Ua3 
In  terms  of  the  velocity  ratio  one  obtains: 


Pa 

K 


1  +  YiMi3[l 


Ha 


/Ua  \a 
ui 


"I 

J 


where,  from  reference  (12): 


v  7R  +  2Ev'(Ti) 

1  5R  +  2E;'(r!) 


where 


(23) 
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It  may  be  seen  in  equation  (23)  that  for  our  range  of  Mach 
number  (10  <  M  <  20)  giving  values  of  temperature  greater  than 
50°R  and  less  than  250or  the  value  of  Ev'  is  zero  and  Yx  -  1.4. 

One  also  needs  the  equation  of  state  and  the  enthalpy 
equation : 


P  -  pRT 


(7) 


H 


7 

2 


RT  +  Ev 


(3) 


The  iteration  procedure  is  carried  out  in  the  following 

way : 


One  first  assumes  ua/ux  -  0  and  calculates  Pa  and  Ha  from 
equations  (21)  and  (22).  From  a  plot  of  equation  (3)  (fig.  11), 
ra  is  obtained  and  pa  calculated  from  (5).  Because  of  the 
continuity  equation,  a  new  value  of  ua/ux  is  given  by  pi/pa. 

The  procedure  is  then  repeated  using  the  new  value  of  ua/ux  until 
two  successive  values  of  this  ratio  differ  by  less  than  the 
desired  accuracy.  Usually  three  or  four  iterations  are  sufficient. 
This  final  value  of  ua/ux  is  used  in  equation  (21)  to  obtain 
Pa/Px. 


c.  In  equation  (20),  the  ratio  poa  /Pa  is  obtained  in  a 

similar  manner  to  the  way  in  which  Pox/Pi  was  obtained  in  the 
equilibrium  nozzle  flow  calculations.  That  is,  ouce  Pa  and  Ta 
are  known  from  the  iteration  procedure,  the  value  of  S  is  calcu¬ 
lated  from  equation  (9).  This  value  is  kept  constant  (isentropic 
flow)  and  the  pressure  at  the  stagnation  point  P^  is  found  from 
equation  (9)  using  this  constant  value  of  S. 

Equation  (20)  is  then  finally  used  to  obtain  the  stagnation 
pressure  ratio  Poj/Pox  •  This  ratio  has  been  plotted  for  each 
case  and  various  stagnation  conditions  versus  Mach  number 
(figs.  12,  13,  and  14)  and  versus  area  ratio  (fig.  15).  In 
figure  14,  all  the  results  for  the  range  3960  *  T0  *  5400  and 
80  atm  *  P o  *  100  atm  fell  on  the  two  curves  shown. 


DISCUSSIOM 

The  purpose  of  these  calculations  as  stated  in  the  intro¬ 
duction  was  to  provide  some  information  on  the  effect  of 
vibrational  energy  on  the  nozzle  flow  in  NOL  Hypersonic  Tunnel 
No.  4.  It  was  hoped  that  this  information  would  assist  in  the 
interpretation  of  experimental  results.  Basically,  the  problem 
is  this:  once  the  stagnation  pressure  ratio  Po*/Pox  ls  experi¬ 
mentally  determined,  how  can  one  obtain  the  static  pressure, 
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static  temperature,  free  stream  velocity,  and  density  in  the 
hypersonic  nozzle  with  vibrational  relaxation  present? 

In  attempting  to  interpret  the  results  found  from  our 
calculations,  one  encounters  the  complication  that  even  in  an 
lnvlscld,  one-dimensional  flow  with  relaxation  present  the  usual 
quantities,  P,  p,  T,  and  u  are  not  uniquely  related  to  the  Mach 
number.  See,  for  example,  figures  12,  13,  and  14.  However,  in 
figure  15  one  sees  that  the  parameter  A/A*  has  succeeded  in 
collapsing  any  differences  between  the  frozen  flow  and  the 
partial  equilibrium  cases  to  a  single  curve.  Hence,  it  seems 
appropriate  to  discuss  the  results  in  terms  of  area  ratio.  But 
it  must  be  noted  that  in  a  physical  flow  there  is  always  a 
boundary  layer  and  one  never  really  knows  what  the  area  ratio  is. 
Furthermore,  a  physical  flow  in  a  conical  nozzle  is  not  one- 
dimensional.  Thus,  we  shall  think  of  the  area  ratio  as  an 
"effective"  area  ratio,  with  no  real  physical  significance,  but 
which  serves  as  a  convenient  quantity  to  plot  quantities  against. 
Its  real  utility  is  as  an  intermediate  parameter. 

Vith  this  in  mind,  we  may  give  an  answer  to  the  question 
posed  in  the  first  paragraph.  Ve  know  that  case  2  is  closest  to 
reality,  but  it  requires  solution  of  the  rate  equation  for  the 
vibrational  energy.  On  the  other  hand,  figure  14  tells  us  that 
when  plotted  versus  area  ratio,  the  ratio  of  pitot  to  supply 
pressure  is  the  same  in  the  completely  frozen  flow  case  as  in 
relaxing  flow.  It  is  of  interest,  therefore,  to  estimate  how 
much  error  is  introduced  when  the  pressure,  temperature,  and 
velocity  are  calculated  using  the  "frozen  flow  area  ratio." 

This  area  ratio  is  obtained  from  figures  like  15,  or  equivalently, 
from  frozen  flow  tables,  such  as  reference  (10).  As  far  as  the 
pressure  is  concerned,  figure  3  shows  that  the  error  is 
negligible.  For  temperature,  from  figure  4,  it  can  be  as  much 
as  9  percent,  whereas  figure  5  shows  that  the  error  can  be  as 
much  as  5  percent  for  velocity.  It  should  be  noted  that  these 
figures  are  given  on  the  basis  of  only  a  limited  number  of  cal¬ 
culations.  It  is  noted,  furthermore,  from  figures  12,  13,  and 
14,  that  for  Mach  numbers  around  15,  the  error  introduced  in 
Mach  number  calculations  using  the  frozen  flow  assumption  is 
between  1  and  2  percent. 


CONCLUSIONS 

In  order  to  determine  the  effects  of  vibrational  energy  in 
the  modified  NOL  Hypersonic  Tunnel  No.  4,  ve  have  considered  the 
flow  to  be  first  in  equilibrium,  then  relaxing,  and  finally 
completely  frozen.  Plots  of  pressure,  temperature,  and  velocity 
versus  area  ratio  for  each  of  these  three  cases  show  that  in 
going  from  frozen  to  equilibrium  flow,  these  three  variables 
increase  at  any  particular  station  in  the  nozzle.  If  one 
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considers  the  relaxing  case,  it  is  found  that  the  flow  remains  in 
equilibrium  up  to  the  throat,  then  freezes  shortly  downstream, 
remaining  frozen  throughout. 

The  flow  around  a  pitot  probe  has  been  examined  for  these 
three  cases,  and  it  is  found  that  when  the  ratio  of  pitot  to 
supply  pressure  is  plotted  versus  area  ratio,  the  relaxing  flow 
case  gives  identical  results  to  the  frozen  flow  case.  Therefore, 
the  authors  recommend  that  the  frozen  flow  area  ratio  be  used  to 
determine  the  flow  variables  in  the  nozzle  once  the  stagnation 
pressure  ratio  is  known.  If  this  is  done,  for  the  cases 
considered,  the  percentage  error  in  assuming  frozen  rather  than 
relaxing  flow  is  estimated  to  be  negligible  for  pressure,  a 
maximum  of  9  percent  for  temperature,  and  5  percent  for  velocity. 
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FIG.  I  VIBRATIONAL  ENERGY  AS  A  PERCENTAGE  OF  INTERNAL  ENERGY  AND  ENTHALPY 
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FIG.  2  COMPRESSIBILITY  FACTOR  FOR  NITROGEN 
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FIG.  4  STATIC  TEMPERATURE  VERSUS  NOZZLE  AREA  RATIO 
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FIG  5  FREE  STREAM  VELOCITY  VERSUS  NOZZLE  AREA  RATIO 
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FIG.  6  VIBRATIONAL  RELAXATION  TIME  FOR  NITROGEN 
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FIG.  8  FREEZING  TEMPERATURE  AS  PERCENTAGE  OF  SUPPLY  TEMPERATURE 
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FI6  9  NUMBER  OF  COLLISIONS  REQUIRED  FOR  NITROGEN  VIBRATIONAL  EQUILIBRIUM 
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FIG  10  TRANSIT  TIME  COMPARED  TO  RELAXATION  TIME  IN  SHOCK  LAYER 
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FIG.  12  STAGNATION  PRESSURE  RATIO  VERSUS  FREE  STREAM  MACH  NO 
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FIG  13  STAGNATION  PRESSURE  RATIO  VERSUS  FREE  STREAM 
MACH  NO. 
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FIG.  14  STAGNATION  PRESSURE  RATIO  VERSUS  FREE  STREAM  MACH  NO 


FIG  15  STAGNATION  PRESSURE  RATIO  VERSUS  NOZZLE  AREA  RATIO 


